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ABSTRACT 

We present a study of the M83 cluster population, covering the disc of the galaxy 
between radii of 0.45 and 4.5 kpc. We aim to probe the properties of the cluster 
population as a function of distance from the galactic centre. We observe a net decline 
in cluster formation efficiency (F, i.e. the amount of star formation happening in 
bound clusters) from about 26% in the inner region to 8% in the outer part of the 
galaxy. The recovered F values within different regions of M 83 follow the same F 
versus star formation rate density relation observed for entire galaxies. We also probe 
the initial cluster mass function (ICMF) as a function of galactocentric distance. We 
observe a signihcant steepening of the ICMF in the outer regions (from —1.90 ±0.11 
to —2.70 ± 0.14) and for the whole galactic cluster population (slope of —2.18 ± 0.07) 
of M 83. We show that this change of slope reflects a more fundamental change of the 
’truncation mass’ at the high-mass end of the distribution. This can be modelled as a 
Schechter function of slope —2 with an exponential cut-off mass (Me) that decreases 
significantly from the inner to the outer regions (from 4.00 to 0.25 x 10^ Mq) while the 
galactic Me is Ri 1.60 x 10^ Mq. The trends in F and ICMF are consistent with the 
observed radial decrease of the S(H 2 ), hence in gas pressure. As gas pressure declines 
cluster formation becomes less efficient. We conclude that the host galaxy environment 
appears to regulate 1) the fraction of stars locked in clusters; 2) the upper mass limit 
of the ICMF, consistently described by a near-universal slope —2 truncated at the 
high-mass end. 

Key words: galaxies:spiral - galaxies:star clusters - galaxies:star formation - 
star: format ion 


1 INTRODUCTION 

Young star clusters are potentially bright, long lived, tracers 
of star formation within their host galaxies. It is therefore 
of paramount importance to u nderstand what govern s the 
cluster formation process (e.g., iLongmore et al.|l2014h and 
how the host galactic e nvironment may affect th eir forma¬ 
tion and evolution (e.g.. lAdamo fc Bastian Il2015l ). 

Statistically we can describe the cluster formation 
as a stochastic process. Numerous works report of an 
initial cluster mass and luminosity function (ICMF and 
ICLF, respectively) being described with good approxima¬ 
tion by a single power-law slope of index close to —2 (e.g., 


Biket^ 200^ de Grhs et al.l l2003l : IWhitmore et al.l l2007l . 
20141 : IChandar et al.ll2014l ). The stochastic nature of clus¬ 


ter formation combined with the power-law distributions of 
the cluster mass generate the so called size-of-sample ef¬ 
fects. These effects are well described by the positive r ela- 
tions between the number of clusters JWhitmorel 1200(3) or 
the youngest brightest cluster in V band, versus 

the g alactic star formation rate (SFR: Earse nl l2002l : feastianl 
I2OO8I) . In other words, galaxies with higher SFR have more 
numerous cluster populations, therefore they are more likely 
to sample the ICMF at higher mass (luminosity) ranges. 
The size-of-sample effect also applies to the time-binning 
of the cluster population. Massive clusters are statistically 
more likely to be older because the galaxy has had a longer 
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J ieriod of time to s ample the high mass bins of the ICMF 
Hunter et alJlioOSl i. 

Increasing evidence, however, suggests that this is not 
the full picture. A steepening at the bright luminosity bins 
of the ICLF has wide ly been observed and reported (e.g. , 
Whitmore_et_^ 19991: Larsen 200^ 5ieles et ^ l2()06lj : 


Bastian et al.ll2012l : IWhitmore et al.ll20l3i 7 This steepening 


could be explained by the pr esence of a truncati on at the 
high mass end of the ICMF dCieles et al.l I 2 OO 6 IJI . Indeed, 
the presence of a truncation mass may explain why we do 
not observe very massive you ng star clust ers close to 10® 
Mq mass in the Milky Way dLarsenlbood i. However, it is 
difficult to establis h the presence of a truncation due to low 
number statistics (iLarse nl l200dl . Several works report an 
ICMF compatible with a Schechter function of slope —2, 
and an exponential cutof f at masses above a certain char¬ 
acteristic mass, Me (e.g ., I^rsen|[2009l : iBastian et al.|[^012l : 
iKonstantopoulos et aklboish . Iiiterestinglv it has been ob¬ 
served by these authors that Me may change as function of 
the host galactic environment. Such an environmental de¬ 
pendence would be expected theoretically due to the varia¬ 
tion of the maximum mass scale for gravitational instability 
in galaxy discs (the "Toomre mass”, Toomre 1964), which 
should lead to an increase of the maximum cluster mass with 
the gas pressure (Kruijssen 2014). 

More evidence of the role played by the galactic en¬ 
vironment on the properties of their cluster population has 
indirectl y been suggested by the so called Tl(U) versus SFR 
relation ijLarserJl20n2l : lAdamo et al.ll201lh . Tl(U) is the frac¬ 
tion of U band light locked in star clusters with respect to the 
total U band light of the galaxy. This fraction is observed 
to increase as function of the SFR of the galaxy, suggest¬ 
ing that the fraction of star formation happening in clusters 
(her eafter F or cluster format ion efhciency, CFE) is increas¬ 
ing l| Adamo fc Bastian 11201,^ . Indeed direct evidence of a 
varying CFE as function of the average galactic star for¬ 
mation rate densi ty (Ssfr) has been observed in numer¬ 
ous galaxies (e.g., |G pddMrLgtaDl2Qlij; lAdamo et al ][20i3; 
lAnnibali et al.ll201ll : Rvon et al.||2014|r am ong m any others). 

In the model proposed by KruiissenI (l2012lf a variation 
in the CFE is directly linked to the star and cluster for¬ 
mation process itself. This model describes star and clus¬ 
ter formation continuously across the density spectrum of 
the hierarchically structured ISM, in which the highest star 
formation efficiencies (and thus bound stellar fractions) are 
reached in the highest density peaksQ Integration of these 
local physics over the gas density spectrum then results in 
a CFE that increases with the gas surface density or pres¬ 
sure, as well as the SFR density (Kruijssen 2012). Given 
the density gradient observed in disc galaxies, the model 
thus predicts that the CFE decreases with the galactocen- 
tric radius. Indeed, using a small fract i on of the whole cluster 
population of M 83, ISilva-Villa et al.l ll2013l . hereafter SV13) 
report evidence of such a variation of the CFE as function 
of galactocentric distances. 

The aim of this article is to probe the imprints that the 


^ In addition, tidal heating by encounters with giant molecular 
clouds (GMCs) imposes a minimum density for the formation of 
bound clusters, which in practice is only important in the very 
highest-density environments. 


galactic environment leaves on the properties of its cluster 
population. The results will shed light on how the star clus¬ 
ter population is a natural outcome of the physics driving 
galactic star formation, and will also help to use star clus¬ 
ter populations as tracers of the SFH of their host galaxies. 
To achieve this goal we study the cluster population of the 
nearby (~ 4.5 Mpe, corresponding to a distance modulus 
of 28.28 mag) spiral galaxy M 83. This galaxy hosts a rich 
cluster po pulation which has been widely analysed in re¬ 
cent years dChandar et ^ 201^ SiJvaWilIa&Jjarsen[ 2011; 


BasHanet_aJj_201l|;_ WIiitnmreet_^ 2011 |_|Bastianet_^|y 


201^^SjjyaWffiajl^l.1 l2013l . l2014l : ichandar et al.l l2014l : 

Hollvhead et al. I l2015h . Several studies suggest that the cur¬ 


rent ongoing starb urst at the centre of this galaxy (e.g., 
IWofford et al.ll201ll') was mostly lik ely produced by a recent 
minor merger i Knaoen et ahlEoiOl ). Outside this region the 
most active star-forming regions in the galaxy are at the 
end of the bar (see Figure [IJ. However, we notice that star 
formation is occurring over the entire spiral arm system and 
in spurs extended in side the inter-arm regions, as in M 51 
llChandar et al]l201lll . 

Only recently M83 has been fully observed with the 
Hubble Space Telescope (HST) with broad and narrow imag¬ 
ing bands covering fro m the UV to the NIR and several 
optical emission lines dBlair et al.ll20l3 'l. Previous studies 
of the cluster population have theref ore focused on the two 
point ings FI and F2 (see Figure [T]). IChandar et al.l (l201Cll . 
l2014h report of a similar ICMF and cluster disruption rate 
in both inner (FI) and outer (F2) fields. They notice a steep¬ 
ening of the ICMF slope and a milder disruption rate in F2, 
but they conclude that within the uncertainties the results 
support a cluster formation and evolution scenario indepen¬ 
dent of the galactic environment. Wit h a similar clus t er cat ¬ 
alogue (comparison p re sented in both iBastian e t al.l ll2012ll : 
Ichandar et al.l (I2014 H. [Bastian et al.l ( 20 111 . 20121 . hereafter 
B12) point out the differences in the photometric properties, 
ICMF and ICLF, disruption rate of the cluster populations 
in the two pointings. The differences cannot be accounted for 
by a recent, significant change in the SFH, so are therefore 
most likely linked to the cluster formation and disruption 
properties of the local galactic environment. Such a change 
in cluster properties is also supported by the change in the 
CFE as function of distance from the centre of the galaxy 
observed by SV13. 

The access to a full cluster catalogue for M83 now 
allows us to verify and put stronger constraints on pre- 
vious findings. Using the whole M83 cluster catalogue, 
ISilva-Villa et al.] (|2014| . hereafter SV14) find a clear environ¬ 
mental dependency on the strength of the cluster disruption 
with higher cluster disruption in the inner regions of the 
galaxy, as already sugg ested by B12. In a recent publication, 
iHollvhead et al.l (I 2 OI 5 I I has studied the Ho morphology of 
the very young star clusters in M83, reporting that clusters 
are already partially or fully exposed within a few Myr af¬ 
ter their formation. Ryon et al (2015) has focused on the 
analysis of M 83 cluster sizes, important for the dynamical 
evolution of these systems. In the present work we will ex¬ 
pand the ICMF analysis published by B12 and the analysis 
of the CFE by SV13 to the whole cluster population within 
the M 83 disk. All together, these works are presenting an 
unprecedentedly detailed and complete study of the cluster 
population of a nearby galaxy. 


© 2015 RAS, MNRAS 000, [T]-?? 












































































































Cluster formation & environment 3 


The paper is organised as follows. Section 2 describes 
the photometric analysis and the properties of the final clus¬ 
ter catalogue. We also include a description of the ancillary 
data used in the analysis. In Section 3, we present our re¬ 
sults. Section 4 contains a discussion of the results within 
the wider context of cluster and star formation in galaxies. 
A summary of the paper is presented in the Conclusions 
(Section 5). 

2 DESCRIPTION OF THE DATA 

2.1 The cluster catalogue 

This analysis is based on the complete M 83 cluster catalogue 
recently published in SV14. The final cluster catalogue has 
been obtained from two multi-band HST imaging datasets, 
(GO 11360, PI: O’Connell, and GO 12513, PI: Blair). The 
two combined datasets cover within 7 pointings the face- 
on M83 disk up to a distance of 4.5 kpc (3.4') from the 
centre of the galaxy (see Figure [T] and Table [!}. Our anal¬ 
ysis is based on WFC3 imaging in 5 optical bands, F336W 
(hereafter U band), F438W {B band), F555W or F547M (V 
band), F657N (narrow filter centred on Hu), and F814W 
(hereafter I band). Cluster detection has been performed on 
the V band with optimised para meter inputs of the sourc e 
extraction algorithm SExtractor llBertin fc ArnorrMll996ll . 
For each extracted source, aperture photometry has been 
performed in all the bands with an aperture radius of 5 px 
(0.2”4.4 pc) and a sky annulus located at 8 px (0.32”) 
and 2 px (0.08”) wide. Isolated compact clusters, visually se¬ 
lected, were used to estimate and add the amount of missed 
flux lost because of the small aperture size used. Moreover 
small differences in calibration between different pointings 
were corrected as explained in SV14. The final catalogue re¬ 
tains only sources with a concentration index larger than 
1.25 mag and positive detections in 4 filters (UBVI). We 
consider detection in Ho? as additional but not required 
constraint. The Vega zero points used to flux calibrate the 
photometry have been retrieved from the WFC3 instrument 
webpagtQ- Mi lky Way foreground extincti on correction has 
been applied llSchlaflv fc FinkbeineiluOllll . As described in 
SV14, the initial catalogue has been visually inspected and 
each source has been assigned a flag describing the morphol¬ 
ogy of the cluster. Class 1 objects are clusters with sym¬ 
metric and compact light distribution. Class 1 objects are 
most likely bound clusters as tested in B12. Class 2 objects 
show asymmetries in the light distribution, multiple peaks, 
potentially contaminant neighbours. These objects closely 
resemble OB associat ions, therefore there is a fair ch ance of 
them being unbound llCieles fc Porteeies Zwai^l201ll l. Class 

3 objects are spurious detections (bad pixels mostly located 
at the end of the chip, background and foreground sources) 
and have been removed from the final catalogue. Any cluster 
candidate located within a radius of 0.45 kpc has also been 
excluded from the catalogue. We excluded this region be¬ 
cause of the high luminosity gradient which causes a change 
in the completeness of the cluster detection and because of 
the different SFH that this region has experienced (B12). 


^ http://www.stsci.edu/hst/wfc3/phot_zp_lbn 


The number of cluster candidates detected in each field is 
listed in Table [T] 


2.2 Physical cluster properties 

The observed spectral energy distribution (SED) of each ob¬ 
ject, reconstructed from the observed integrated fluxes in all 
the bands, has been compared to stellar evolutionary models 
including treatme nt of nebular gas emiss ion and continuum 
(Yggdrasil models [Zackrisson_et_^]j20^) . The fitting algo¬ 
rithm is described in lAdamo et al.l ll2010ll. The model fluxes 
are reddened accordingly to the ICardelli et all (Il989ll pre¬ 
scription. The algorithm simultaneously fits age and extinc¬ 
tion of the cluster, while the mass is the normalisation factor 
to scale the models to the observed SEDs. Uncertainties on 
the derived ages, masses, and extin ction are estimated ac¬ 
cordingly to the recipe presented in lAdamo et al.l ll2012ll . 

Although final cluster catalogues for FI and F2 were 
available (B12) we decided to repeat the analysis for these 
two fields to make the dataset as homogeneous as possible. 
The comparison with the B12 catalogue shows some minor 
differences in the ages of the clusters between the latter cat¬ 
alogue and the new one as reported in SV14. The difference 
was produced by an erroneous extinction value systemati¬ 
cally applied to the I band which was a factor of 1.26 higher 
than the correct one. The cluster catalogue, published in 
this work, contains the correct and most updated cluster 
photometry and properties for M 83. As already discussed in 
SV14 and in the present work, the main conclusions reached 
in B12 and SV13 are not affected by this error. 


2.3 Estimates of the star formation rates 

In this work, we have estimated current (between 1 and 10 
Myr) and average (between 10 and 50 Myr) star formation 
rates (cSFR and (SFR), respectively) in different regions of 
the galaxy as follows. The SFR is compared to the cluster 
formation rate (CFR) to derive F. Since clusters have been 
detected in optical wavelengths, we derive SFRs from opti¬ 
cal based tracers. Thus, we use Ha as tracer of the cSFR 
and stellar counts to derive the recent SFH (up to 50 Myr). 
This choice may not take into account the hidden SFR bet¬ 
ter captured by an IR tracer, however, our SFR estimates 
will suffer similar extinction patterns as the clusters do, can¬ 
celling out the effects. 

To estimate the cSFR in different regions of the galaxy 
we used the public ly available Ha continuum subtracted 
image of M 83 (see iDale et ^ l2009l . for details on the re¬ 
duction process and calibration). We calibrated the total 
luminosity estimated from this frame (after masking resid¬ 
uals at the location of foreground stars) to th e total Ha 
luminosity published by iKennicutt et al.l (120081') , giving us 
a correction factor that takes into account [Nil] contamina¬ 
tion and foreground galactic extinction. An average visual 
extinction of Av ~ 0.2 mag has been applied to the total 
Ha flux. To convert Ha luminosity into SFR we use the 
relation by IKennicutt fc EvansI ll2012ll . 

To estimate the (SFR) we derive the recent SFH 
using colour-magnitude diagrams of the resolved stellar 
population. The technique and method are described in 
ISilW-Villa fc LarsenI (I 2 OI 2 I I and SV13. Stellar photometry 
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Figure 1. Left: M83 image in the red visual band {R band), the footprints of the seven HST pointings necessary to cover the bulk 
of the galactic body are overlaid. Centre and right: the Ha continuum subtracted image of MSS shows the regions where the current 
star formation is occurring in the galaxy. Annuli of equal area (centre) and containing equal number of clusters (right) are overplotted. 
The central region, contained inside the green circle , has b een excluded from the analysis. Orientation and resolution of the images are 
indicated on each frame. Image credit at I Dale et al.l <l2003l 'l. taken from NASA/IPAC Extragalactic Database (NED). 


and analysis will be discussed in a forth-coming paper (Silva- 
Villa et al in prep.). For consistency with the cluster and 
Ha analyses, we corrected the final stellar photometry by 
an average visual extinction of 0.2 mag. Because of different 
detection limits between the arm and inter-arm regions we 
applied a conservative approach. Completeness tests show 
that our SFH is robust between 10 and 50 Myr, above and 
below which the significancy of our results drops. The (SFR) 
is thus derived averaging the SFH over 10 to 50 Myr. 

2.4 Estimates of the gas surface density 

To estimate the gas surface density, E(H 2 ), we use the 
velocity-integrated CO (l-O) intensity map published by 
iLundgren et al.l ll2004al ). to which we refer for a descrip¬ 
tion of the dataset and reduction steps. CO intensity was 
derived within each bin used in our analysis. To convert the 
CO intensity to a molecular gas mass, we adopt a conver¬ 
sion factor of Xco = 2.3 x 10^° K km s“^ cm^0 The derived 
values are in g ood agreement with av eraged E(H 2 ) reported 
in their work jLundgren et al.ll2004al ). 


3 RESULTS 

3.1 Analysis approach 

To understand how the galactic environment influences star 
and cluster formation we analyse and compare the cluster 
properties in different regions of M83. The galactic region 
contained between the inner radius of 0.45 kpc and outer 
radius of 4.5 kpc has been divided in four bins. As shown in 

^ We adopt this value for consistency with iLundgren et all 
((SoOl^). This choice is consistent with the Galactic conversion 
factor of Xc.n = 2.0 x 10^^ K km s“^ cm^ with an uncertainty 
of about 30% llBolatto et al.l[2013h . 


the central and right panels of Figure [T] we choose two bin¬ 
ning approaches, equal surface area bins and bins containing 
the same number of objects (the latter approach already ap¬ 
plied in SV13). For each bin we derive the following quanti¬ 
ties (see Table [2]): cSFR (1-10 Myr); (SFR) (10 - 50 Myr); 
r(l - 10 Myr) and r(10 - 50 Myr); predicted CFE (F**'); 
gas surface density, E(H 2 ); luminosity of the youngest most 
massive cluster, mass of the most massive cluster (in 

the age range 3 to 100 Myr), Mmax', characteristic mass, Mr,; 
and slope of the mass function if a single power-law is used, 
Q. In Figure[21 we show the age-mass diagram of the cluster 
population contained in each bin of equal area. We include 
both class 1 (red dots) and 2 (black dots) objects. We do not 
see any significant difference within the two classes and each 
bin has similar fractions of classl versus class 2 objects. The 
blue and green boxes show the cluster population used to 
derived the CFE as described below. The age distributions 
of Figure [3 clearly show the size-of-sample effect, i.e. the in¬ 
creasing Mmax as a function of age (and also of decreasing 
Mmax as a function of galactocenric distance). 

The CFE is defined as the ratio between the CFR and 
the SFR within the same age interval. We derive the CFR 
in the same interval used for the SFR (1-10 and 10-50 Myr). 
The CFR is the ratio between the total stellar mass in clus¬ 
ters and the corresponding age range. The current CFR is 
estimated from clusters of age between 1-10 Myr and clus¬ 
ter mass above a thousand Mq. This quantity is compared 
to the cSFR estimated within the same bin to derive r(l-10 
Myr). Similarly, the average CFR is estimated for clusters 
with age between 10 and 50 Myr and mass above 5000 Mq 
and directly compared to the (SFR) thus deriving F (10-50 
Myr). To derive the total CFR we assume an underlying 
ICMF power-law of slope ?2, with lower and upper mass lim¬ 
its of 100 Mq and two times the maximum observed cluster 
mass in the region under investigation, respectively. The to¬ 
tal CFR will be the result of the observed stellar mass in 


© 2015 RAS, MNRAS 000, [T]-?? 























Cluster formation & environment 5 


Table 1. Description of the HST dataset used in this analysis. All the data have been taken with the WFC3. The exposure time of each 
band is included in brackets. N^i is the final number of clusters (class 1 & 2) in each field. We include within brackets how many of these 
clusters belong to class 1. 


Field 

PI 

Filters 

N.i" 

FI 

O’Connell 

F336W (1890 s), F438W (1880 s) 
F555W (1230 s), F657N (1484 s) 
F814W (1213 s) 

1387 (672) 

F2 

O’Connell 

F336W (2560 s), F438W (1800 s) 
F547M (1203 s), F657N (1484 s) 
F814W (1213 s) 

1098 (547) 

F3 

Blair 

F336W (2579 s), F438W (1799 s) 
F547M (2682 s), F657N (1799 s) 
F814W (1379 s) 

1142 (595) 

F4 

Blair 

F336W (2589 s), F438W (1809 s) 
F547M (2682 s), F657N (1809 s) 
F814W (1389 s) 

1243 (687) 

F5 

Blair 

F336W (2589 s), F438W (1809 s) 
F547M (2682 s), F657N (1809 s) 
F814W (1389 s) 

1199 (676) 

F6 

Blair 

F336W (2579 s), F438W (1799 s) 
F547M (2682 s), F657N (1799 s) 
F814W (1379 s) 

1414 (774) 

F7 

Blair 

F336W (2579 s), F438W (1799 s) 
F547M (2682 s), F657N (1799 s) 
F814W (1379 s) 

859 (521) 


® The clusters located within 0.45 kpc from the centre have not been in¬ 
cluded in this final catalogue; 


clusters more massive than the observational imposed limit 
and the missing fraction. 

The two mass limits applied at different age ranges to 
estimate the CFR are well within the detection limits of 
our sample. They are chosen to be substantially higher than 
the detection limit to avoid completeness issues. However, 
these mass limits may cause possible misclassifications due 
to stochastic effects in the stellar initial mass function (IMF) 
as pointed out by several studie s (e.g., IPopescu fc HansonI 
I2OI0I : l^uesneau fc Lanconll2010h . In a upcoming work, our 
deterministic approach to derive cluster properties is com¬ 
pared to the predictions obtained by SLUG (Krumholz et al 
in prep). SLUG is a bayesian fitting code t hat uses stochas¬ 
tic IM F models to derive cluster properties (iKrumholz et al.l 
I2OI5I I. Krumholz et al (in prep) find significant misclassifica- 
tion of the cluster mass below 1000 M©, i.e., the determinis¬ 
tic approach used in this work significantly underestimates 
cluster masses below a 1000 Mq. 

We estimate the fraction of missing stellar mass in clus¬ 
ters with mass below the two thresholds, assuming a power- 
law s lope distribution of —2 and a lower mass limit of 100 
Mq (iLongmore et al.ll2014l L Only class 1 clusters are used 
to derive CFRs. Due to this conservative approach and the 
possible effects of stochasticity, our GFE determinations can 
be considered lower limits to the real values. Monte Garlo 
simulations of cluster populations are used to estimate the 
uncertainties associated with the GFR (and hence F) val¬ 
ues. We take into account a 0.1 dex error in both derived 
age and mass of each cluster (uncertainties from the SED 


fitting technique). We also include the poissonian likelihood 
of recovering the observed GFR. This uncertainty accounts 
for the observed number of clusters within the considered 
mass and age range as randomly drawn from a power-law 
cluster mass function with upper mass limit equal to two 
times the observed maximum mass value in each bin. 

In Table [ 2 ] one can see that the SFR and F derived with 
different methods agree within 3 sigma in each bin, confirm¬ 
ing that the SFR, as well as the GFE, has not changed sig¬ 
nificantly in the last 50 Myr (SV14). A different binning of 
the data does not produce any major impact in the derived 
quantities. As already discussed in SV13, bins containing the 
same number of clusters should mitigate the size-of-sample 
effect. In the M83 case, differences between equal surface 
area and equal number of cluster bins are mostly restricted 
to the regions containing a very active site of star-formation 
for the M83 disk, i.e. the edges of the bar (see the central 
and right panels of Figure [T]). Therefore, even if local vari¬ 
ations between the two binning approaches are observed in 
the central bins covering the end of the bar region, the radial 
trends discussed in the next sections are preserved. 


3.2 On the fraction of star formation happening 
in clusters as function of distance within the 
galaxy 

In Figure [3 we show how F changes with galactocentric dis¬ 
tances for bins of equal area. The decreasing trend observed 
here confirms the finding reported by SV13, but for a sub- 
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Table 2. Galactic and cluster population properties derived for the whole population, in annuli of equal area, or containing equal number 
of clusters. Only class 1 clusters have been used in this analysis. The typical uncertainty on is a factor of two (see Section 3.2). 

The Me, ot have been derived for classl clusters younger than 100 Myr and more massive than 5000 Mq. 


annulus 

kpc 

cSFR 

Moyr-i 

(SFR) 

M0yr-l 

pi — lOAfyr 

% 

plO—SOAfyr 

% 

pth 

% 

E(H2) 

M0pc ^ 

Mbright 

mag 

Mmax 

lO^M© 

Me 

10® M0 

a 

galactic 

0.45-4.50=^ 

0.82±0.8 

1.20±0.12 

18.2±3.0 

15.Oil.9 

19.0'^ 

47.7 

-11.60 

1 qk+0-45 
i.ya_Q 29 

1.60i0.30 

-2.18i0.07 

equal area 

0.45-2.30 

0.25±0.02 

0.21±0.02 

19.3±4.0 

26.5i4.0 

28.6 

62.4 

-11.60 

1 

J--4:0_o.i9 

4.00i0.80 

-I.90i0.ll 

2.30-3.20 

0.31±0.03 

0.43±0.04 

26.0±3.5 

19.2i2.6 

20.4 

34.9 

-11.47 

l.95t°;« 

1.00i0.20 

-2.20i0.ll 

3.20-3.90 

0.15±0.02 

0.34±0.03 

13.1±2.5 

9.8il.6 

13.5 

19.4 

-10.68 

0.59t°;« 

0.55i0.11 

-2.20i0.12 

3.90-4.50 

O.lliO.Ol 

0.22±0.02 

13.7±3.1 

8.0il.5 

12.6 

17.7 

-10.55 

0.29t°;°^ 

0.25i0.05 

-2.70i0.14 

equal number 

0.45-2.30 

0.26±0.03 

0.25±0.02 

18.9±4.4 

23.2i4.2 

28.6 

62.4 

-11.60 

1 

J--4:0_o.i9 

SO.OOi 10.00 

-I.90i0.ll 

2.30-2.60 

0.12±0.01 

0.25±0.02 

34.0±5.1 

20.8i4.0 

26.3 

52.6 

-11.47 

l.95t°;« 

1.50i0.30 

-2.05i0.12 

2.60-3.30 

0.21±0.02 

0.37±0.04 

19.8±3.6 

10.8i2.1 

17.4 

27.6 

-10.27 

0.59t°;« 

0.60i0.12 

-2.35i 0.13 

3.30-4.50 

0.23±0.02 

0.33±0.03 

14.0±2.6 

12.6i2.2 

13.0 

18.3 

-10.55 

0.391°;°^ 

0.30i0.06 

-2.50i0.13 


^ Properties of the galaxy and cluster population integrated within 0.45 and 4.5 kpc. 
^ Value from iKruii^e 3 1 I 2012 I') . 



Figure 2. Cluster age versus mass diagram of each bin of equal 
area. The black dashed lines show the detection limits in the V 
band {My ^ —5.5 mag in bin 1 and My ^ —5.3 mag in the 
other 3 bins), the filter used to extract the source catalogue (see 
main text). Black dots are cluster candidates of class 2, red dots 
clusters belonging to class 1. The blue and green delimited areas 
show the clusters used to derive the current and averaged CFR. 
The limit on the cluster mass is well above the detection limits 
and has been imposed to mitigate the effect of stochasticity of 
the IMF. 

stantially larger cluster sample. The two T values (for each 
age range) follow each other quite closely. F decreases in¬ 
side out from about 26% to 8%. On different scales (right 
y-axis), we overplot the gas surface density, S(H 2 ) esti¬ 
mated within the same regions. E(H 2 ) shows a similar de¬ 
cline (black dashed line) suggesting a link between the CFE 
and the amount of molecular gas available for star forma¬ 
tion. The plot includes al so the prediction s for F obtained 
with the fiducial model bv iKruiissenI (I2ni2h . The model pre¬ 
dicts a higher absolute value than the observed one in the 
centre but agrees quite well with the values we find in the 


disc, within the uncertainties. We notice that the i n ner p oint 
(green star symbol), taken from iGoddard et al.l (120101) . is 
likely to be underestimated due to incompleteness caused by 
the strong differential extinction in this region. In the disc 
(between 0.45 and 4.5 kpc), the observed relative decline of 
F is closely reproduced by the model. 

The good agreement between the predicted and ob¬ 
served gradients of F as a function of galactocentric distance 
suggests that the Kruijssen cluster formation model can be 
used to understand the cluster formation process. The model 
predicts that the formation of bound stellar clusters takes 
place in the highest-density peaks of a hierarchically struc¬ 
tured ISM. At these high densities, the gas goes through 
a large number of free-fall times on a short time-scale and 
therefore reaches a high star formation efficiency. As a re¬ 
sult, the expulsion of any residu al gas by feedback doe s not 
unbind the stellar distribution dKruiissen et al.l l2012l) and 
a bound cluster survives. At lower densities, the low star 
formation efficiencies yield unbound associations^ Putting 
these local considerations in a galactic context, the model 
predicts that clusters form most efficiently at high gas pres¬ 
sures (and hence gas surface densities), because these condi¬ 
tions lead to higher density peaks and thus favour bound 
cluster formation. In summary, the model thus predicts 
that the CFE is set by galactic-scale properties such as 
the gas surface density and, therefore, the star formation 
rate density Ssfr throug h the Schmidt-Kennicutt relation 
llKennicutt &: Evansll2012 ). 

The two panels in Figure 4 support the above inter¬ 
pretation. We show the CFE as a function of Esfr for our 


An additional included effect at high gas surface densities 
(S > 10^ Mq pc“^) is that tidal perturbations by dense gi¬ 
ant molecular clouds unb ind lower-density clus ters that otherwise 
would haye suryived (cf. iKruiissen et al.ll201lh . 
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Figure 3. Cluster formation efficiency as function of galacto- 
centric distance. F has been estimated for two age ranges and 
different mass limits as indicated in the inset. Only clusters clas¬ 
sified as class 1 have been used for the analysis. Therefore, the 
CFEs can be considered lower limits to the real fractions. Ha 
emission is the tracer used to estimate the SFR between 1 and 10 
Myr, while direct stellar counts have been used to derive the SFR 
in the age range 10-50 Myr. The horizontal bars show the width 
of each annulus. The area within each annulus is the same. The 
green star shows the F value for the star burst region confined i n 
the centre of the galaxy and reported bv iGoddard et ^ ] | |2010D . 
The right-hand y-axis shows the molecular gas surface density 
(S(H 2 )). See text for more details. 

spatially resolved data points in M 83 as well as for a compi¬ 
lation of literature measurements of the CFE in entire galax¬ 
ies. Both the entire galaxies and the different regions within 
galaxies follow the same trend of increasing CFE with the 
SFR density. We include the CFEs estimated for the two 
types of M 83 binning (equal area and equal cluster number 
bins). We also estimate the CFE in each of the 7 fields shown 
in Figure 1. We show the M 83 CFEs for both age ranges 
(left and right panel). Overall, the agreement with Kruijs- 
sen’s fiducial model (black dotted line) is very encouraging 
and adds to the evidence supporting a tight link between 
galactic environment and stellar cluster formation. 

3.3 Cluster mass function as function of distance 
within the galaxy 

B12 reports evidence of a steepening of the ICMF in the 
cluster population of F2 relative to that of FI (see the po¬ 
sition of the two fields in the galaxy in Figure [T]). The two 
cluster populations are better fitted by a Schechter mass 
function of slope —2 and Mc~ 1.5 x 10® and 0.5 x 10® Mq 
in the inner and outer fields, respectively. Access to a com¬ 
plete cluster catalogue for this galaxy enables us to further 
investigate how the galactic environment possibly influences 
the shape of the ICMF. Using the same analysis technique 
as in the previous section, we construct cumulative distri¬ 
butions of the cluster mass of each radial bin. Cumulative 
distributions are favoured to avoid any subjectivity in the 


choice of the mass bin size and because they are more sen¬ 
sitive to variations. Statistically, this method has the same 
significanc e of distributions of bins contain ing same number 
of objects llMafz Apellaniz fc Ubedall2005h . with our distri¬ 
butions containing only one object per bin. In the previous 
paragraph, to estimate F we have limited our analysis to 
cluster ages overlapping the age range to which our chosen 
SFR traces are sensitive. However, to build the ICMF we 
do not need to apply the same limits. Similarly to the limits 
used by B12 and in other recent works, we apply a lower 
mass limit of 5000 Mq and an age cut between 3 and 100 
Myr. For comparison, we generate Monte Carlo simulations 
of cluster populations with the same cluster numbers of the 
observations. We assume the same lower mass cutoff as in 
the observations (5000 Mq) and an upper limit to the clus¬ 
ter mass of 10^ Mq, that is much larger than the observed 
maximum cluster mass so it approximates the condition of 
no upper mass truncation of the ICMF. 

We discuss below the results obtained using only class 1 
clusters. In Figure ET] of the appendix we include the same 
analysis but building the ICMFs of both class 1 and 2 ob¬ 
jects, instead. The recovered absolute values vary within a 
factor of 2. Therefore, the observed trends and conclusions 
remain unchanged confirming that our results are not driven 
by the selection criteria. 

In FigureO we show the observed and simulated ICMFs 
for the entire cluster population (top row) and for the clus¬ 
ters contained in bins of equal area (second top row to bot¬ 
tom). The results do not change if cluster populations in 
bins of equal numbers are used instead (see Table El). Each 
panel shows the observed cumulative distribution (red dots), 
the median (black solid line) and the limits containing 50 
(dashed lines) and 90 % (dotted lines) of the 2000 cumula¬ 
tive ICMF realisations. We use different functions, i.e., single 
power law of slope —2 and no upper mass limit; single power 
law of some best-fitting slope that is not necessarily —2; a 
Schechter function of power-law slope —2 and exponential 
truncation above a variable characteristic mass, M^. 

To verify the goodness of the agreement between the ob¬ 
served and simulated distributions we use the Kolmogorov- 
Smirnov statistic. The resulting probability, p(KS), that the 
two distributions are actually produced by the same parent 
population are included in each panel. In the first column 
we attempt to reproduce these distributions with a single 
power-law function of slope —2. The innermost bin cluster 
mass distribution (top left panel) is consistent with a power- 
law function of slope —2. The significance of the deviations 
from this function increases in the outer bins and the whole 
M83 cluster population (panels in the first column). The 
middle column panels show an attempt to reproduce the 
ICMFs with single power laws of steeper slopes. We derive 
the value of the slope which maximises p(KS). Although the 
upper mass distributions are not well reproduced by a single 
power-law function, the KS test favours steeper slopes than 
the initial value —2. 

The change in the slope from ~ —1.9 to —2.7 from 
the inner to the outermost bin is quite striking. In Ta¬ 
ble El we report the uncertainties associated with the slope 
in each radial bin as well as for the whole galaxy. We de¬ 
rive the uncertainties on the slopes determining the index 
interval which contains 68% of the Monte Carlo simula¬ 
tions for each cluster sample. As widely discussed in pre- 
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log ^SFR [Moyr"' kpc"^] 
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Figure 4. Fversus SFR surface density for class 1 clusters. The same diagram shows the position of the cluster formation efficiency 
estimated within different regions of the galaxy (see top inset) using circular bins of same area, annuli containing the same number of 
clusters, the position of each HST pointing (7 fields). On the left plot, T has been derived from clusters younger than 10 Myr, on the 
right plot, from clusters in the age range 10 to 50 Myr. Literature data and predictions have been included for comparison (see Table IBTI 
in the Appendix for references and a complete list of the values used in this plot). 


vious works (e.g., iGieles et al.ll2006lj) . a steepening in the 
upper mass (luminosity) end of a distribution is often com¬ 
patible with th e presence of a tr uncation. In the third col¬ 
umn we use a ISchecht^ (Il976l ) function, i.e. dN/dM oc 
M‘^ exp { — M/Mc), with slope a — —2 and variable trun¬ 
cation mass Me, as indicated in the inset of each plot (see 
also listed values in Table [2]). When using a Schechter func¬ 
tion rather than a single power law slope —2, the agreement 
between simulations and observations improves for bin 2, 
bin 3, bin 4 and for the cluster population altogether. This 
shows that in 80% of the cases, a Schechter function gives 
a substantially better match to the observed ICMF than a 
single power-law fit. A unique feature of the Schechter func¬ 
tion is the ability to reproduce the observed distributions at 
the high-mass ends. Statistically a steeper power-law func¬ 
tion and a Schechter type function give similar solutions, 
although the latter seems to be preferred everywhere but 
the innermost bin. Only there, the KS statistic favours a 
single power-law £t with a slope slightly less steep than —2, 
even though this slope cannot be ruled out within the uncer¬ 
tainties. For this reason, the Schechter £t produced with a 
power-law slope of index —2 does not improve the KS prob¬ 
ability. If instead of using a Schechter function with slope 
—2 we use —1.9 the KS probability becomes p(KS)= 0.58, 
which is still lower than the p(KS)=0.794 obtained for the 
single power-law fit. In section 3.4, we discuss the physics 
driving this difference. 


3.4 The most massive cluster as function of 
distance within the galaxy 

Widely used in the cluster literature, the brightest clus¬ 
ter luminosity versus the SFR of the gala xy is a rela¬ 
tion dominat ed by the size-of-sample effect llLarse nl l2002l: 
Bastianll2008l). T his relation, as well as the one presented by 
Whitmore! ( 2^000| ) . shows that cluster formation is a stochas¬ 
tic process, i.e., the mass or luminosity of the most massive 
cluster in a galaxy increases with the SFR, because a larger 
number of formed clusters in the population increases the 
probability of forming more massive or luminous clusters. 
In Figure|6]we show the SFR versus relation includ¬ 

ing all the data available in the literature. We also add the 
brightest young cluster recovered in each of the 4 bins of 
M83 (pur ple dots ). The dashed black line is a fit to the 
sampl e by iLarse 3 1I2002I) and presented in IWeidner et al.l 
(I2OO4I) . Using Monte Carlo simulations under the assump¬ 
tion that only a relative small f raction o f star formation is 
happening in clusters (r~ 8 %). iBastianI (I2OO8I ) reproduces 
the observed trend. If all star formation would happen in 
clusters then the observed points should follow the dotted 
line on the left side of the dashed one. The scatter of the data 
points is quite large. Many factors can affect the location of 
cluster samples in this diagram ^ such as a recent change in 
the SFH of the galaxyJlBastim gOO^, or a difference in the 
F of each galaxy llAdamo et al l201lj) . 


Our results suggest that the ICMF changes signihcantly 
as function of distance from the centre and that the ICMF 
can be overall described by a power-law function of slope 
—2 and a decreasing Me as function of distance from the 
centre of the galaxy. This finding is in agreement with the 
analysis presented by B12 for a smaller cluster sample in 
M83. It clearly suggests that the change in the star forma¬ 
tion properties of the galaxy as function of radius is affecting 
the shape of the upper mass-end of the ICMF. 


Even though the SFR versus relation is mainly 

driven by the size-of-sample effect, i.e. the stochastic sam¬ 
pling of the most massive and luminous clusters from the 
ICMF, this does not rule out that the galactic environment 
still plays an important role in shaping the properties of the 
cluster population (see Section 4). In Figure (5] we include 
new Monte Carlo simulations of cluster populations formed 
with a Schechter ICMF of different Me assuming a fixed 
cluster formation efficiency of 8% (red solid lines). The pres¬ 
ence of a truncation mass in the ICMF clearly introduces a 
bend in the simulated SFR versus relation. However, 
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Figure 5. Cumulative ICMFs (red filled dots) of the whole population and of the clusters within bins of same area (from top to bottom). 
Only clusters classified as 1 (see Figure 1 in the Appendix for the same plots but class 1 Sz 2 instead), more massive than 5000 Mq 
and ages between 3 and 100 Myr have been included. Median (solid line), quartiles (dashed line) and extended boundaries (dotted line) 
of Monte Carlo simulations of cluster populations with the same number of objects as the observed ones are included in each panel. 
In the first column, we assume a pure power-law ICMF with slope —2. In the middle column we still assume a power-law ICMF but 
with varying slopes (see insets), chosen to maximise the Kolmogov-Smirnov probability that the observed and simulated distributions 
are drawn from the same parent distribution, p(KS). In the right column we assume that the ICMF is better described by a Schechter 
function of slope —2 and a exponential cut-off. Me, chosen to maximise p(KS). 
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log(SFR) [Me/yr] 


Figure 6. The luminosity of the brightest cluster plotted against 
the star-formation rate (SFR) of the host galaxy. This plot con¬ 
tains a compilation of all data available in the literature. See Ta¬ 
ble in the Appendix for a detailed description of the datasets. 
The dash ed line is the b est fit to the sample of galaxies plotted as 
triangles llLarse i] |2002t) . The dwarf galaxy s ample is represented 
by green filled stars fsee lAdamo Sz Bastian 1120151 . for a descrip- 
tion of the sample). Squares are the sample added by iBastianI 
ll2008l') . who derived the expected SFR vs. M^^^^^relation if 100% 
of stars are born in clusters with a power-law ICMF with index 
—2 (black dotted line). Blue stars are the sample of luminous blue 
compact galaxies studied bv lAdamo et al] ll201ll) . The green hor¬ 
izontal bar sho ws the range of S FR of dwarfs which do not have 
young clusters llCook et al 1[2 oT 1 . The M 83 brightest young clus¬ 
ter of each the four bins of equal area are plott ed as purple cir¬ 
cles. T he orange diamonds show the sample from lWhitmore et alj 
||2014| ~). The red solid lines show Monte Carlo simulations of clus¬ 
ter populations formed with a Schechter CMF with different Me 
and a cluster formation efficiency of 8%. 

this bend is absent in the observed SFR versus rela¬ 

tion, because galaxies with higher SFRs tend to have higher 
truncation masses Me (see below). 

The presence of a varying upper mass cutoff in the 
ICMF is supported by o bservations of GMCs i n local galax¬ 
ies. In the M 51 galaxy, IColombo et al.l (l2014l l report clear 
differences in the number, mass, luminosity, density, and 
mass distributions of GMCs within different environments of 
the galaxy (e.g. inter-arm and arm regions, molecular ring, 
etc.). The mass of the most massive CMC in M51 changes 
as fun ction of the ga lactic environment. With a similar argu¬ 
ment, [harse^ ll2009ll noticed that the maximum GMC mass 
observed in the Antennae system is larger than the GMG 
mass found in local spiral galaxies. Clusters form within 



Figure 7. The mass of the most massive cluster found in each bin 
of same area is plotted as function of galactocentric distance (red 
triangles). The shadowed area in orange show the uncertainties 
on the mass estimates. The Me derived in Figure [5] are plotted 
using cyan dots, while the blue band show 20% tolerance margin 
for the derived values. The p redicted cluster My are included as 
black stars (see the text and lKruiiss^l2014h . 


GMCs, and their final mass is a fraction of the total mass of 
the GMC. The total fraction of gas mass converted into stars 
is usually described by the star formation efficiency, which 
observations of galactic star-forming regions suggest is onl y 
of a few per cents (e.g.. lElmegreen|[2002l : lEvans et al.ll2009l ). 
The maximum fraction of the GMC mass that can end up 
in a single, most massive cluster then requires thi s star for¬ 
matio n efficiency to be multiplied by the CFE llKruiissenI 
l2014h . The observation of an ICMF truncation mass at the 
high-mass end is possibly linked to the limit imposed by 
the galactic environment to the maximum mass that GMCs 
can have. Theory predicts that this mass-scale should corre¬ 
spond to the maximum mass able to overcome the galactic 
differ ential rot ation and unde rgo collapse (i.e. the iToomrel 
1 19641 mass, see lKruiissenl[20l3 and references therein). 

In Figure [7] we show the observed maximum cluster 
mass within each bin in M83 (red triangles and mass un¬ 
certainties derived from SED fitting), the recovered Me in¬ 
cluding 20% uncertainties, and pre dicted ma x imum cluster 
mass-scale using equation 3 from iKruiissenI ll2014h (black 
symbols). It is important to stress here that although we 
refer to Me as a truncation it is not a sharp limit. The trun¬ 
cation of a Schechter function reflects an exponential decay 
of the probability distribution at the high-mass end. The 
chances of having clusters more massive that Me are lower 
than if a pure p ower-law fun ction of slope —2 are used but 
not impossible dLarse nl l2009l ). Indeed, we see that the ob¬ 
served M-max and the recovered Me are in agreement within 
a factor of 2. To derive the predicted maximum cluster mass 
we use the averaged F (mean of the values at 1-10 and 10- 
50 Myr ) observed in each bin of equal area and we as¬ 
sume that the GMG mass changes between 2 and 0.8 x 10^ 
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M(70 (iLundgren et alJ [2004bl ’) in each bin and fix the star 
formation efficiency to 5%. In spite of the approximations, 
the predicted and the observed trend agree remarkably well, 
supporting the role of the galactic environment in shaping 
the upper mass limits of the ICMF. 


4 DISCUSSION 

In this work, we probe the impact of the environment on the 
cluster formation process. Our analysis reinforces previous 
findings and presents unambiguous imprints left by the M 83 
galactic environment on its own cluster population. 

Cluster formation is a stochastic process, however the 
efficiency of the process, i.e. the amount of star formation 
locked in bound star clusters, varies as function of the Esfr. 
Several works have shown this link on g alactic scales (e.g., 
iGoddard et al.ir20ld : lAdamo et al.l[201l[ among many oth¬ 
ers). This effect is also observed on sub galactic scales, when 
looking on differ ent regions of the same galaxy (e.g., SV13, 
iRvon et al.ll2014l i. and our new analysis confirms this result. 
We find that the cluster formation process is tightly linked 
to the gas surface density, just li ke the SFR density is via 
the S chmidt-Kennicutt relation llSchmidtlll959l : iKennicrittl 
Il998h . As the gas pressure increases (and hence the gas or 
SFR surface density) the density spectrum of the ISM broad¬ 
ens and a larger fraction of the density fluctuations reaches 
densities high enough to collapse on short time-scales and 
achieve high star formation efficiencies. This locally more ef¬ 
ficient star formation enables the young stellar st ructure to 
rema in bound when any residual gas is expelled dKruiisserJ 
l2012lb As a result, the CFE increases with the gas pressure, 
gas surface density, and through the Schmidt-Kennicutt re¬ 
lation also with the SFR density. 

However, F appears to not grow indefinitely as it obvi¬ 
ously cannot attain values exceeding 100%. At gas densities 
above E(H 2 ) > 1000 Mq/pc”^ (which corresponds to r~70 
%) F will not increas e significantly (see dotted line in Fig¬ 
ure [T] |Krui]sse3[2Ql2)- Denser gas favours star and cluster 
formation but at the same time destroys low density systems 
more easily via CMC encounters. These tidal perturbations 
imply that a CFE of 100% is never reached. 

The presence or absence of a truncation at the high- 
mass end of the ICMF has been long deba ted in t he lit - 
erature. With simple numerical simulations iLarse 3 (l2006li 
shows why it is so challenging to prove the presence of a 
truncations. For example, in Section we find that the 
ICMF of the innermost bin can be well described by a sin¬ 
gle power-law mass function of slope —1.90 ±0.11. However, 
the presence of a truncation cannot be ruled out from the KS 
test. Large cluster populations are necessary to well sample 
the high mass end of the ICMF, but in local galaxi es clus- 
ter populations have usually a few hundred objects llLarserJ 


® We prefer using the observed maximum cloud masses rather 
than calcula ting the Toomre mass u sing the velocity dispersion 
profiles from iLundgren et al.l ll2004bl i. These authors discuss that 
beam smearing yields velocity dispersions which (especially at 
small galactocentric radii of < 2 kpc) can be overestimated by 
a factor of 2. This would subsequently lead to severely overesti¬ 
mated Toomre masses (Mx oc cr^) by up to an order of magnitude. 


l2006h . Some previous works in the past have found evi¬ 
dence of a possible truncation a nd/or variation of the ICMF 
within the same galaxy fe.g.. iLarse 3 I 2 OO 9 I : iBastian et al.l 
I 2 OI 2 I : iKonstantopoulos et al.ll2013ll or steepening of the c lus- 

ter luminosity function f recently I Whitmore et aLll2014l L In 

our work we observe a clear steepening of the ICMF, well 
represented by a Schechter type distribution of slope —2, 
and a variable Me as function of distance from the centre of 
the galaxy. The Schechter function combines the power-law 
slope and the presence of a limit above which the proba¬ 
bility to form more massive clusters goes exponentially to 
zero. It is not a sharp limit. Interestingly, we observe a clear 
decline of Me within the same galaxy. As for the CFE, the 
decrease of the truncation mass is most likely related to the 
gradient of the gas pressure (and hence gas and SFR surface 
density). The galactic environment thus poses a direct limit 
on the maximum possible mass-scales of GMCs and stellar 
clusters. 

The globular cluster mass function has been reported 
to be well fitted by an evolved Schechter function by 
Ijordan et al.l (l200itl . This function takes into account the 
stellar mass loss suffered by globular clusters at the low-mass 
end of the mass function as well as the steep, non-gaussian 
decline at the high-mass end. The Me of each globular clus¬ 
ter population increases with the total B band luminosity 
(mass) of the host galaxy. Dynamical friction alone is not 
able to reproduce the observed decline of Me for lower lu¬ 
minosity galaxies suggesting a link between Me and the host 
galaxy at the moment the globular cluster mass function was 
created. 

The results reported in our analysis of M 83 and sup¬ 
ported by evidence found in previous works can shed light 
on the cluster formation process not only in local galax¬ 
ies but also at high redshift. The young star cluster mass 
function can be described by an almost universal ICMF of 
slope —2 and a variable Me which appears to be a func¬ 
tion of the host galaxy environment, in that it increases 
with the gas pressure. Combining this with the aforemen¬ 
tioned increase of the truncation mass of the globular clus¬ 
ter mass function with the host galaxy mass, we infer that 
the gas pressure at the time of globular cluster formation 
likely increased with the host galaxy mass , as is indeed pro¬ 
posed by recent models (^Kruiissenll2015^ . Our results can 
help to understand the physical process under which glob¬ 
ular clusters formed. Both young and ancient cluster pop¬ 
ulations can potentially be used as tracers of the star for¬ 
mation process of their host galaxies. However, it is impor¬ 
tant to keep in mind that the physical conditions of the 
ISM where globular clusters have formed are quite different 
and nowadays observed only in merging starburst systems 
(e.g., iHoltzman et al.l [l992l : I Whitmore et al.l 1 19931 1. More¬ 
over globular cluster populations have most likely survived 
because during merging events they have been relocated in 
the halos of their galaxies, thereby esca ping the gas-rich , dis¬ 
ruptive bodies of their host galaxies (lKruiissenll2015li . For 
instance, the M83 cluster population will be affected by clus¬ 
ter disruption according to the results of SV14. In the ab¬ 
sence of major events which could move clusters out of the 
disk, this population will not survive for a Hubble time. 
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5 CONCLUSIONS 

The aim of this work is to probe the link between cluster for¬ 
mation, the observed statistical cluster properties and the 
galactic environment where they form and interact. With 
respect to previous analyses, we now have access to a com¬ 
plete cluster catalogue covering the vast majority of the M 83 
galaxy. We have estimated the cluster formation efficiency 
across the whole disc of M83 using several techniques and 
SFR tracers. We have sliced the cluster catalogue in 4 bins 
either containing the same number of clusters (as done by 
SV13) or having the same area. Bins containing the same 
number of clusters remove the size-of-sample effect. How¬ 
ever, we show in Section |3] that our results are not affected 
by the binning technique. 

In general, we find that, within a factor of two, the SFR 
derived using Ha as tracer or stellar count techniques are in 
very good agreement. Therefore, we conclude that the M 83 
star formation rate has been nearly constant in at least the 
last several tens of Myr as suggested in B12, SV13, and 
SV14. 

We derive lower limits to the cluster formation effi¬ 
ciency, r as function of galactocentric distance and using 
the position of the 7 pointings across the galactic body. The 
CFE of the whole cluster population of M83 has a lower 
limit of ~ 18 %. With a wider coverage of the galaxy and 
thus a more numerous cluster population we are able to put 
strong constraints on the derived F values. We see a net 
decline in F from about 26% in the inner bin to 8% in the 
outer bin. A similar decline is observed in the averaged gas 
surface densi ty and in the pr edictions made with the fidu¬ 
cial model bv iKruiissenI (l2012l l. We see that the derived val¬ 
ues of r versus Esfr in different regions of M 83 follow the 
F-Esfr relation of galaxies. We conclude that the relation 
appears to hold not only on galactic but also on sub-galactic 
scales. The dependence of F on Esfr arises from a more fun¬ 
damental dependence on the gas pressure (or surface den¬ 
sity) through the Schmidt-Kennicutt relation. Therefore, the 
amount of stars locked in clusters appears to be regulated by 
the same mechanism that regulates the overall star forma¬ 
tion process. Similarly to the Schmidt-Kennicutt relation, a 
well-calibrated T-Esfr or r-E(H2) relation can potentially 
be used to make realistic predi ctions of cluster formation 
in cosmological simulations le.g. lKruiissenl[2012l l. However, 
more effort is needed it to empirically link cluster properties 
to star formation and GMC properties. 

Another aspect that can potentially reveal imprints of 
the galactic environment on the cluster population is the 
value of the possible truncation in the upper mass-end of the 
ICMF. To investigate this issue we have built the ICMFs of 
the cluster population contained in each galactocentric ra¬ 
dius bin. We observe a significant steepening of the ICMF 
as function of galactocentric distances. Monte Carlo simula¬ 
tions of a single power-law function without upper mass limit 
are not able to reproduce the observed ICMFs. The steep¬ 
ening can only be consistently reproduced If a Schechter 
function of slope —2 and exponential cut-off at the high- 
mass end with a varying Me is used. Me decreases signif¬ 
icantly in the outer radial bins. This finding is consistent 
with the observed decrease of the E(H 2 ) and likely of the 
maximum CMC mass. The probability that the ICMF of 
M 83 is drawn by a single power-law mass function of slope 


—2 is ~ 4 X 10“^. The probability increases significantly if 
a steeper power-law function {a « —2.2) is used. The best 
representation of the M 83 ICMF is a Schechter function of 
slope —2 and Mc« 1.5 x 10® M©. Overall we conclude that 
the young star cluster mass function can be described by an 
almost universal ICMF of slope —2 and a variable Me which 
appears to be a function of the host galaxy environment. Up¬ 
coming systematic surveys of stellar cluster populations will 
enabl e us to quantify an d constrain this variation further 
(e.g., ICalzetti et al.l[2015ll . 
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APPENDIX A: THE CMF OF CLASSl & 2 
CLUSTERS IN M83 

We report here the CMF analysis of classl and 2 objects 
contained in bins of equal area. 
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Figure Al. Cumulative ICMFs (red filled dots) of the whole population and of the clusters within bins of same area (from top to 
bottom). Both classl and class 2 systems more massive than 5000 Mq and ages between 3 and 100 Myr have been included. Median 
(solid line), quartiles (dashed line) and extended boundaries (dotted line) of Monte Carlo simulations of cluster populations with the same 
number of objects as the observed ones are included in each panel. In the first column, to generate the fake populations we assume a pure 
power-law ICMF with slope —2. In the middle column we still assume a power-law ICMF but with varying slopes (see insets), chosen to 
maximise the Kolmogov-Smirnov probability that the observed and simulated distributions are drawn from the same parent distribution, 
p(KS). In the right column we assume that the ICMF is better described by a Schechter function of slope —2 and a exponential cut-off, 
Me, chosen to maximise the KS probability, as indicated in the inset. 
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Table Bl. a) lAdamo et alJ (blue star symbols in Figure [bl; b) iLarse 3 ll20n2l~l (black triangles in Figure 0 , only a fraction of 

the original catalogue is listed while the values of the omitted ones are listed with the most recent estimated values available in the 
literature; c) The values used for each binning (purple filled dots in Figure [6)1 are reported in Table[2]of the main text; d) iBasti^ l|2008l') 
(black squares in Figure!^: ellWhitmore et all ll2014li (orange diamonds in Figure O, some of the galaxies in the original list have been 
omitted because iLarserjFoO^ presented data with a larger coverag e of the galaxy than in t his paper. The dwarf gala xies (green star 
symbo ls i n Figur^bjl include da tapo ints from the following pap ers: rijArinibajj_et a3_j |2niih j_QlJPasgimli et^l.l ll2nilh: fe llOoddard et ^ 
ll2010l l ; i) lAnnibali et alJ ll2009ll ; j ) iRafelski fc Zari tskvl : fcl | Johnsoir'er'ahril200Cll h 1) ICook etan ~ il2012l l: ml lLim fc Led ll2015l h ^ 

In Figure |4] we show the combined F derived bv' ICook'er'aLf 1 I 2 OI 2 I 1 using all the clusters younger than 10 Myr detected in their dwar f 
sa mple. The SFR de nsity reported is not a single value but a range enclosing all the single galactic values, o') ISilva- Villa &: LarsenI l(20l3); 
o) IR von et al.l ll2ni4l ~l: r) This work. Values obtained considering the age range 1-10 Myr, see Table 2 in the main text for more details. 


Galaxy 

Mbright 

[mag] 

SFR 

Mgyr-l 

Esfr 

MQyr^^^Kpc”-^ 

F 

% 

ESO 338“ 

-15.50 

3.20 

1.55 

50.0±10.0 

Haro 11“ 

-16.16 

22.00 

2.16 

50.0113° 

ESO 185-1013“ 

-14.36 

6.40 

0.52 

26.0±5.0 

MRK 930“ 

-15.00 

5.34 

0.59 

25.0±10.0 

SBS0335-052E“ 

-14.28 

1.30 

0.95 

49.0±15.0 

NOG 247'’ 

-10.20 

0.0360 

0.18x10“*^ 

- 

NGC OOO*” 

-9.90 

0.0779 

0.49x10-3 

- 

NGC 1156*’ 

-11.10 

0.1842 

3.07x10-3 

- 

NGC 2403** 

-9.90 

0.3376 

0.97x10-3 

- 

NGC 2835** 

-10.90 

0.0920 

0.73x10-3 

- 

NGC 2997*’ 

-12.90 

1.8604 

3.07x10-3 

- 

NGC 2997‘J 

- 

- 

9.4x10-3 

10.0±2.6 

NGC 3184*’ 

-10.60 

0.3956 

1.72x10-3 

- 

NGC 3621*’ 

-11.90 

0.8801 

1.67x10-3 

- 

NGC 5204*’ 

-9.60 

0.0315 

0.83x10-3 

- 

M83*’’“ 

-11.70 

2.2842 

13.76x10-3 

- 

M83 (centre)*’ 

- 

0.39 

0.54 

26.7ll'3 

M83 (outer)P 

- 

0.39 

0.013 

5.6±0.6 

M83 (0.45-4.5 kpc)’’ 

-11.60 

0.82 

0.013 

18.2±3.0 

NGC 5585*’ 

-10.80 

0.0336 

0.32x10-3 

- 

NGC 6744*’ 

-11.00 

0.4309 

0.62x10-3 

- 

NGC 6946*’ 

-13.00 

2.5392 

4.60x10-3 

- 

NGC 6946'’(WFPC2) 

- 

0.17 

4.60x10-3 

12.511:1 

NGC 7424*’ 

-11.40 

0.1728 

0.18x10-3 

- 

NGC 1741*’ 

-15.00 

4.8819 

12.78x10-3 

- 

NGC 5253*’ 

-11.10 

0.2114 

7.29x10-3 

- 

IC 1613*’ 

-5.80 

0.0004 

0.05x10-3 

- 

LMC 

-11.U 

0.1201*’ 

1.52x10-3*’ 

5.8±0.5*’ 

NGC 4214*’ 

-12.04 

0.0798 

3.80x10-3 

- 

DDO 50*’ 

-7.91 

0.0108 

1.26x10-3 

- 

DD0168*’ 

-7.58 

0.0043 

0.85x10-3 

- 

DDO 165*’ 

-8.34 

0.0005 

0.18x10-3 

- 

Sextans 

-7.12 

0.0037 

2.29x10-3 

- 

NGC 3521*’ 

-11.50 

1.3676 

3.58x10-3 

- 

NGC 4258*’ 

-12.60 

0.9926 

0.70x10-3 

- 

NGC 5055*’ 

-11.40 

1.5019 

2.98x10-3 

- 

M51*’ 

-12.80 

4.7454 

8.21x10-3 

- 

NGC 7252^* 

-13.40 

5.40 

- 

- 

NGC 6240’* 

-16.40 

140.00 

- 

- 

NGC 2207’* 

-13.60 

2.20 

- 

- 

NGC 1275’* 

-15.30 

12.40 

- 

- 


APPENDIX B: A COMPLETE COMPILATION OF THE DATA SAMPLES AVAILABLE IN THE 
LITERATURE 

Complete compilation of the data available in the literature and plotted in Figure and [4] The footnote associated with the 
name of the galaxy gives the reference for all the values reported on the corresponding row, unless a different note on the 
values specify a different source. Some galaxies are reported more than once, with different values of SFR and Esfr, because 
r and related quantities (SFR, and Ssfr) have been derived using only partial coverages of galaxies. 
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Table Bl. - Continued 


Galaxy 

Mbright 

SFR 

SsFR 

F 


[mag] 

Meyr-l 

MQyr-^Kpc--^ 

% 

M82'^ 

-14.80 

7.00 

- 

- 

NGC 3597^* 

-13.30 

10.80 

- 

- 

IRAS 19115'^ 

-16.80 

192.00 

- 

- 

NGC 1533'^ (Al) 

-7.17 

0.37x10-3 

- 

- 

NGC 1533'^ (A2) 

-5.71 

0.25x10-3 

- 

- 

NGC 1533'^ (A3) 

-6.16 

0.18x10-3 

- 

- 

NGC 2623*^ 

-14.50 

51.00 

- 

- 

NGC 3256^* 

-15.70 

46.00 

0.62 

22.91^:1*' 

NGC 7673*^ 

-14.70 

4.90 

- 

- 

NGC 6745"^ 

-15.00 

12.20 

- 

- 

NGC 1140'* 

-14.80 

0.80 

- 

- 

Milky Way*' 

- 

0.1508 

0.012 


NGC 45= 

-10.83 

0.12 

- 

- 

NGC 45P 

- 

0.05 

1.02x10-3 

5.2±0.3 

NGC 406= 

-11.75 

0.29 

- 

- 

NGC 628= 

-11.84 

0.23 

- 

- 

NGC 1300= (FI) 

-11.00 

0.30 

- 

- 

NGC 1300= (F2) 

-11.52 

0.32 

- 

- 

NGC 1309= 

-13.80 

1.70 

- 

- 

NGC 1313= 

-10.98 

0.22 

- 

- 

NGC 1313P 

- 

0.68 

0.011 

3.2±0.2 

NGC 1483= 

-10.01 

0.11 

- 

- 

NGC 3627= 

-11.97 

0.37 

- 

- 

NGC 4038= 

-15.25 

2.43 

- 

- 

NGC 4394= 

-10.25 

0.14 

- 

- 

NGC 4395= 

-9.79 

0.07 

- 

- 

NGC4395P 

- 

0.17 

4.66x10-3 

l.OitO.6 

NGC 4736= 

-10.44 

0.04 

- 

- 

NGC 5055= 

-9.61 

0.02 

- 

- 

M101= (FI) 

-11.38 

0.25 

- 

- 

M101= (F2) 

-11.57 

0.20 

- 

- 

NGC 6503= 

-10.51 

0.07 

- 

- 

NGC 7793= 

-9.65 

0.07 

- 

- 

NGC 7793P 

- 

0.15 

6.51x10-3 

2.5±0.3 

NGC 4449-'' 

-12.60 

1.0 

0.04 

9.0 

NGC 1569 

-13.9*’ 

0.369 

0.03® 

13.9±0.8'' 

NGC 1705' 

-13.8 

0.31 

0.046 

- 

SMC-^^ 

-9.94 

0.043 

0.001 

A 9+0.2/1 

^•^-0.3 

He2-10'= 

-12.5 

0.20 

0.20 

- 

NGC 2366* 

-8.52 

0.094 

2.4x10-3 

- 

UGC4305* 

-8.88 

0.114 

2.2x10-3 

- 

UGC 4459* 

-7.9 

0.003 

1.0x10-3 

- 

UGC 5336* 

-8.59 

0.004 

0.00037 

- 

IC 2574* 

-9.12 

0.11 

1.2x10-3 

- 

UGC 5692* 

-9.44 

7.34x10-3 

0.5x10-3 

- 

UGCa292* 

-7.46 

5.15x10-3 

0.01 

- 

UGC 8201* 

-9.16 

31.6x10-3 

2.1x10-3 - 


UGC 9128* 

-5.45 

0.3x10-3 

0.3x10-3 

- 

UGC 9240* 

-5.90 

5.24x10-3 

1.0x10-3 

- 

Dwarf sample® 

- 

- 

[0.01, 10.0] x 10-3 

5.0 

ICIO"* 

-10.4 

0.07 

0.03 

4.2 
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